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Abstract. The photoproduction of and n°n ± pairs off 40 Ca has been investigated with the TAPS 

detector using the Glasgow photon tagging spectrometer at the Mainz MAMI accelerator. Data have been 
taken for incident photon energies in the energy range from 200 - 820 MeV. Total cross sections have been 
extracted from threshold up to the maximum photon energy and invariant mass distributions of the pion 
pairs have been obtained for incident photon energies between 400 - 500 MeV and 500 - 550 MeV. The 
double 7T° invariant mass distributions show some relative enhancement with respect to the mixed charge 
channel at small invariant masses. The effects are smaller than previously observed for lead nuclei and 
the distributions do not significantly deviate from carbon data. The data are in good agreement with the 
results of recent calculations in the framework of the BUU model, with careful treatment of final state 
interaction effects but without an explicit in-medium modification of scalar, iso-scalar pion pairs. This 
means that for Ca most of the experimentally observed effect can be explained by final state interactions. 
Only at low incident photon energies there is still a small low mass enhancement of the double ty° data 
over the BUU results. 

PACS. I3.60.Le meson production - 25.20.Lj photoproduction reactions 



1 Introduction 

In-medium properties of hadrons are a hotly debated topic 
since they are closely related to the properties of low- 
energy, non-perturbative QCD. QCD at high energies or 
short scales (r < O.lfm) is well explored by perturbative 
methods. However, at larger distances the perturbative 
picture breaks down and at scales between 0.1 fm and 1 
fm the full complexity of the interaction manifests itself in 
the many body structure of hadrons composed of valence 
quarks, sea quarks, and gluons. Chiral symmetry is at the 
very heart of the theory. In the limit of vanishing current 
quark masses the QCD Lagrangian is invariant under chi- 
ral rotations, right and left handedness of quarks is con- 
served and right and left handed fields can be treated inde- 
pendently. The explicit breaking of the symmetry due to 
finite u,d current quark masses (5-15 MeV) is small. How- 
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ever, it is well known that spontaneous breaking occurs 
since the ground state, the QCD vacuum, has only part of 
the symmetry. This is connected to a non-zero expectation 
value of scalar qq pairs in the vacuum, the chiral conden- 
sate. The symmetry breaking is clearly reflected in the 
hadron spectrum. Without it, hadrons should appear as 
mass degenerate parity doublets, which is neither true for 
baryons nor for mesons. The first J' K =l/2~ excited state 
in the baryon spectrum, the Sn(1535) lies much above the 
J 7r =l/2+ nucleon ground state. The lowest lying J*=l~ 
vector meson, the p, has a smaller mass than the J 7T =1 + 
a±, and the J n =0~ pion (the Goldstone boson of chiral 
symmetry), is much lighter than its chiral partner, the 
JTV—Q+ (j-nieson. 

However, model calculations (see e.g. [1]) indicate a 
significant temperature and density dependence of the chi- 
ral condensate which is connected with a partial restora- 
tion of chiral symmetry at high temperature and/or large 
densities. An observable consequence is a density and tem- 
perature dependence of hadron masses. These effects will 
be most pronounced for the high temperature/density con- 
ditions probed in heavy ion reactions, but they should be 
already significant at zero temperature and normal nuclear 
density as probed by pion and photon beams. In- medium 
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modifications of vector mesons have been searched for via 
the spectroscopy of di-lepton pairs emitted in heavy ion 
reactions by the CERES experiment [2[[3] and more re- 
cently the NA60 collaboration [4] , which reported a strong 
broadening of the in-medium spectral function of the p- 
meson. A signature for in-medium modifications of p and 
lo mesons was also reported from the E325 experiment in 
12 GeV p+A reactions at KEK [5J. An in- medium mass 
shift of the lo meson has been found in photon induced 
w-production from heavy nuclei [6] . Due to their coupling 
to vector mesons, nucleon resonances should also be mod- 
ified in the nuclear medium. Self-consistent calculations of 
the respective spectral functions [7] predict in particular 
significant effects for the Di3(1520) resonance, which have 
also been searched for in photon induced meson produc- 
tion reactions [M91 ITUIIT1"] . however so far without conclu- 
sive results. As discussed in detail by Lehr and Mosel [12] . 
part of the problem is, that due to the averaging over the 
nuclear volume, a possible broadening of the resonance 
can be only observed via decay channels which cause the 
broadening, but not in other exclusive reactions. 

A particularly interesting case is the mass-split be- 
tween the J^CP pion and the J 7r =0 + a- meson. The naive 
assumption that the two masses should become degener- 
ate in the chiral limit is supported by model calculations. 
A typical result is the density dependence of the mass 
calculated in the Nambu-Jona-Lasino model by Bernard, 
Meifiner and Zahed [13] . The er-mass drops significantly as 
function of nuclear density. Compared to the vacuum the 
predicted effect is already large for normal nuclear mat- 
ter density po at which the pion mass stays still stable. 
The nature of the a meson itself has been controversely 
discussed in the literature. The review of particle prop- 
erties [13] lists as a the /o(600) with a mass range from 
400 - 1200 MeV and a full width between 600 MeV and 
1000 MeV. Recently, Caprini, Colangelo and Leutwyler 
|15j have derived precise predictions for mass and width 
from dispersion relations. In some approaches it is treated 
as a pure qq (quasi)bound state [T3 l [T6 l fT7] . in other ap- 
proaches as a correlated wir pair in a / = 0, J 71 ' = + state 
[18 l fT9 l l20] . But in any case a strong coupling to scalar, 
iso-scalar pion pairs is predicted. As a consequence, the 
different model approaches agree in so far, as they predict 
a significant in-medium modification of the invariant mass 
distribution of such pion pairs. This is either due to the 
in-medium spectral function of the a meson |16j or the in- 
medium modification of the pion-pion interaction [19] due 
to the coupling to nucleon - hole, A - hole and N* - hole 
states. In any case the predicted net effect is a downward 
shift of the strength in the invariant mass distributions of 
scalar, iso-scalar pion pairs in nuclear matter. 

First experimental evidence for this effect has been 
reported by the CHAOS collaboration from the measure- 
ment of pion induced double pion production reactions 
[2"Tll22[|23[|2"4Tl2"5] . The main finding was a build-up of strength 
with rising mass number at low invariant masses for the 
7r + 7r~ final state, but not for the 7r + 7r + channel where 
the <7 as neutral particle cannot contribute. A similar ef- 
fect was found by the Crystal Ball collaboration at BNL. 



They observed a low-mass enhancement of strength for 
heavy nuclei in the ir~ A — > Air°ir° reaction [26] . 

Pion induced reactions have the disadvantage that the 
signal is diluted by final (FSI) and initial (ISI) state inter- 
action of the pions, so that effectively only the low-density 
surface zone of the nuclei is probed. Pions can be produced 
in the entire volume of the nuclei by photon induced reac- 
tions, although their FSI also suppresses the contributions 
from the deep interior of the nuclei. The FSI can be min- 
imized by the use of low incident photon energies, giving 
rise to low-energy pions which have much larger mean free 
paths than pions that can excite the Z\-resonance [9] . Pho- 
toproduction of the different charge states of pions from 
the free proton and the quasifree neutron has been pre- 
viously studied in detail with the DAPHNE [27 1 [28 l l3o j l3T l 
[32] and TAPS detectors [331I341I351I361I371I38] at MAMI in 
Mainz from threshold to the second resonance region and 
for the 7r°7r° channel also at higher incident photon ener- 
gies at GRAAL in Grenoble [35] (see [3D] for an overview). 



First results from a measurement of double 



and 



7r°7T photoproduction off carbon and lead have been re- 
ported in [41] . A shift of the strength to lower invariant 
masses was found for the heavier nucleus for the ir°ir° 
channel but not for the mixed charge channel. In the 
present paper, experimental details and the results of an 
additional measurement of double pion photoproduction 
off calcium nuclei are presented and compared to model 
calculations. 



2 Experimental setup 



The experimental setup was identical to the measurement 
of double pion photoproduction off carbon and lead nuclei 
reported in [41] . Some further experimental details have 
been discussed in [9]. 




Fig. 1. Setup of the TAPS detector at the Mainz MAMI accel- 
erator. The beam entered the target chamber from the lower 
right edge. 
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Quasi-monochromatic photons were produced with the 
Glasgow tagged photon facility [42] at the Mainz MAMI 
accelerator |43j . The experiment covered the photon en- 
ergy range from 200 - 800 MeV with an energy resolution 
of approximately 2 MeV per tagger channel. The target 
was a 10 mm long cylinder of calcium with 50 mm diame- 
ter and the beam-spot size at the target position was ap- 
proximately 30 mm. The reaction products from the target 
were detected with the electromagnetic calorimeter TAPS 
[44,45]. This detector consists of 510 hexagonally shaped 
BaF2 crystals of 25 cm length with an inner diameter of 
5.9 cm. As shown in fig.[T]they were arranged in six blocks 
of 64 modules and a larger forward wall of 138 modules. 
The blocks were arranged in one plane around the target 
at a distance of 55 cm from the target center and at polar 
angles of ±54°, ±103°, and ±153°, while the forward wall 
was placed 60 cm away from the target center at 0° and 
the photon beam passed through a hole in the center of 
the forward wall. Each detector module is equipped with 
an individual plastic veto detector, read out by a separate 
photomultiplier. The setup covered approximately 37 % 
of the full solid angle. 



3 Data analysis 

The first step of the analysis for both final states - ir°ir° 
and tt ^ - requires the identification of photons. This was 
done with the help of the veto detectors, a time-of- flight 
analysis (time resolution: 550 ps FWHM), and a pulse 
shape analysis of the signals from the BaF2 scintillator. 
The latter is based on the two different components of 
the scintillation light from BaF2 crystal, which have very 
different decay times. The relative intensity of the fast 
component depends on the ionization density produced 
by the incident particle, in particular it is suppressed for 
nuclcons in comparison to photons. This feature was ex- 
ploited by integrating the output signals over a short gate 
(40 ns) and a long gate (2 jus). Both signals E s (short 
gate) and E w (wide gate) were calibrated for photons so 
that in an E s -versus-E w plot photons are lying on the 45° 
line. This is shown in the left hand side of fig. [2] For the 
discrimination between photons and nucleons a represen- 
tation in polar coordinates (see right hand side of fig. [2]) 



is more efficient. The respective coordinates are the ra- 
dius R = ^jE 2 a + ~E\ and the angle <S> = arctan(E s /E w ). 
The combination of the different discrimination methods 
provides a clean sample of photons. 

For the investigation of the 7r 7r° final state events with 
four photons were selected. The invariant masses of all 
possible combinations of the four photons into two pairs 
were calculated and the 'best' combination was chosen via 
a x 2 analysis which minimizes 

with 1 < i / j ^ ( < 4, where is the pion mass, 
the the invariant masses of the photon pairs, 

and the Am la are the associated uncertainties. A two- 
dimensional representation of the invariant masses of these 
pion pairs is shown in fig. [3] Events with both invariant 
masses in the range 100 - 150 MeV were selected. 




Fig. 3. Invariant mass of photon pair 1 versus invariant mass 
of photon pair 2 for events with four detected photons. 




Fig. 2. Pulse shape analysis. Left hand side: plot of E B versus 
E w . Right hand side: plot in polar coordinates. 



At incident photon energies below the r\ production 
threshold this analysis results in a practically background- 
free data sample for the ir°ir° final state. Above the r\ 
threshold background arises from 77 — > 7r n 7r°7r° decays, 
where the decay photons from one ir° escape detection 
due to the limited solid angle of the calorimeter. This 
background can be partly suppressed by a missing mass 
analysis. It is assumed that the reaction occurs quasi-free 
off one participant nucleoli. The initial momentum of this 
bound nucleon is neglected. The mass of the missing par- 
ticle - the participant nucleon - is then determined by the 
energy of the incident photon E^ and the four-vectors of 
the four decay photons apart from the effects of Fermi mo- 
tion, which broaden the missing mass distributions. The 
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Fig. 4. Missing mass spectra for double 7r° photoproduction 
from Ca. Left side: incident photon energy range 400 - 500 
MeV (below rj threshold), middle and right hand side above rj 
threshold: E y around 650 and 800 MeV, respectively. Closed 
circles: experiment, shaded histograms: simulation of 2tt° pho- 
toproduction (around zero) and rj photoproduction, solid his- 
tograms: sum of both. 



missing mass AM is given by: 



AM 
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m N 
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£ 

i=i 



where E~ li and P T 



m N . 



(2) 



are energies and momenta of the decay 
photons and tun is the nucleon mass. The resulting dis- 
tributions are shown in fig.[?]for energies below and above 
the rj production threshold. They are compared to Monte 
Carlo simulations which take into account Fermi motion, 
a simple approximation of final state interaction (FSI) of 
the mesons, and the detector acceptance and detection ef- 
ficiency. The simulation of the detector response was based 
on the GEANT3 package [46] . At incident photon energies 
below the rj threshold the data agree very well with the 
results of the simulation and no indication of background 
is visible. The background from rj decays is prominent 
at higher incident photon energies, but the data can still 
be reproduced by a superposition of simulated double tt° 
and rj — > 7r°7r 7r° events. Due to the large broadening of 
the missing mass structures by Fermi motion and FSI (re- 
sponsible for the shoulders of the double it distributions) 
a clear separation of the two reaction channels is not pos- 
sible. However, the cross section for rj photoproduction is 
simultaneously measured via their rj — > 77 decay. With 
this cross section and the known decay branching ratios 
it is possible to simulate the contribution of 77 — > 7r 7r°7r° 
decays to the analysis of double ir° decays. If this is done 
without any cut on missing mass, the result becomes inde- 
pendent of the simulated shapes of the distributions, but 
the statistical quality suffers since at high incident photon 
energies the rj background contribution is large. Therefore 
as a compromise only events with missing mass between 
-100 and 100 MeV were selected, so that the background 
contribution from the rj decays did not exceed 20 %. 

The analysis of the 7r°7r ± final state is more difficult. 
In the first step events with two photons and the invari- 
ant mass of the photon pair (see fig. [3 left hand side) 
between 100 and 150 MeV were selected. In the second 
step an additional particle was required. Candidates for 
charged pions were then selected by a time-of-flight ver- 
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Fig. 5. Time-of-flight versus energy. The upper part shows 
a three-dimensional representation of the TOF-versus-E spec- 
trum. The bands for electrons, charged pions, protons, and 
even deuterons are clearly visible. The bottom part shows a 
projection. Left hand side: data, right hand side MC simula- 
tion for the 7r°7r ± reaction (only the charged pions were tracked 
by GEANT). The lines indicate the applied cuts. 



sus energy analysis. The corresponding distributions of 
deposited energy E p versus time-of-flight T p are summa- 
rized in fig, [SJ They are compared to the results of an 



MC simulation of 



production. The shape of the 



pion bands in the data and the simulation arc in good 
agreement. It should be noted, that ir + and ir~ behave 
identically in the TOF-versus-E spectra. This is not triv- 
ial since the tt~ get absorbed by nuclei once they have 
deposited their kinetic energy, while the 7r + decay via 
the chain n + — > fJ^v^ — > e + v e u^ after they have been 
stopped. However, the kinetic energy of the decay muons 
is small and their lifetime (2.2 xlO" 6 s) is much longer 
than the electronic gates so that the energy deposition 
from the decay positrons is not measured. In the MC sim- 
ulations particles were only tracked for times correspond- 
ing to the hardware gate lengths used in the experiment. 
Since the time-of-flight is more precisely known than the 
deposited energy of the pions, their kinetic energies were 
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re-calculated from the time-of- flight. The region accepted 
for pion candidates in TOF-versus-E is indicated in fig. [5l 
In order to suppress background from r\ — » 7r°7r + 7r~ decays 
it was furthermore required in the event selection, that no 
second particle fulfilled the charged pion condition. 

The pion band in the TOF-versus-E analysis is not 
background free (see fig. [5]). There is some chance that pro- 
tons leak into it. This happens for example for high energy 
protons which punch through the BaF2 modules and thus 
deposit too little energy or for protons that scatter on their 
way to the detector (e.g. in the target chamber). These 
misidentified protons are in particular a problem at low 
incident photon energies, where the cross section for the 
potential background final state pir° from single ir° pro- 
duction via the A resonance is larger than the ir ^ cross 
section by more than an order of magnitude. A suppres- 
sion of this background is only possible via the reaction 
kinematics. However, one must take into account that pro- 
tons, which are misidentified as pions, are assigned wrong 
kinetic energies since the energies of charged particles are 
based on the time-of-flight measurement using the mass 
of the particle. Therefore, only constraints based on the 
kinematics of the n° mesons can be used. At very forward 
angles this background was so large that the forward wall 
of the detector setup was not used for the charged pion 
analysis, which means that charged pions were only iden- 
tified for laboratory polar angles larger than 22°. 
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Fig. 6. Missing mass spectrum of 7r°7r + events for the back- 
ground hypothesis p-ir (Ej= 400 - 460 MeV). 

The most efficient method for the identification of sin- 
gle 7T° photoproduction is again the missing mass Am 
calculated under the assumption of a quasi-free reaction 
mechanism (analogous to eq.[2]with the sums running only 
over the two decay photons). The result of this analysis 
for the critical low energy region is shown in fig. [6] The 
plot on the left hand side shows the missing mass for the 
hypothesis of single ir° production for events where only 
the detection of a neutral pion was required. It shows a 
clear peak around zero corresponding to single ir° produc- 
tion and the contribution from multiple pion production 
at larger missing masses. The right hand side shows the 
spectrum under the condition that an additional particle 
was detected which fulfills the pion TOF-versus-E condi- 
tion. The data are compared to a MC simulation of the 
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Fig. 7. Left hand side: invariant mass spectrum of photon 
pairs. Right hand side: invariant mass spectrum of 7r°7r ± 7r ± 
events. The dotted line indicates the nominal position of the 
77 peak. The dashed histogram is the mirror image of the left 
hand side of the peak (see text). 



7r°7T final state. Some residual background from misiden- 
tified protons is visible at missing masses lower than 150 
MeV. For further analysis only events with a single ir° 
missing mass above 150 MeV were accepted. 

Finally, the photoproduction of 77 mesons was analyzed 
for the two decay modes rj — > 77 and 77 — » 7T°7T + 7r~. Since 
the cross section for the 40 Ca(7,?y)A reaction is known 
from previous experiments [47], this reaction can be used 
to estimate systematic effects in the detection of photons 
and charged pions and thus it helps to approximate the 
systematic uncertainties of the double pion production 
cross sections. The 77—5-77 decay channel was identified 
with a standard invariant mass analysis of photon pairs 
(see fig. [7J left hand side). Events with invariant masses 
between 500 MeV and 600 MeV were selected. This is the 
range where the invariant mass peak is practically free 
of background, which has been investigated in detail for 
previous 77 production experiments with TAPS at MAMI 
at incident photon energies below 800 MeV ([47 1 148 49 1 
[50] . For the analysis of the 77 — > 7T°7r + 7r~ decay in a first 
step 7t^7t^ pairs were selected with the same cuts as for 
double pion production. In a second step a further hit ful- 
filling the TOF-versus-E condition for charged pions was 
required and finally the invariant mass of the three pions 
was constructed (see fig. [3 right hand side). The spectrum 
is dominated by the invariant mass peak of the 77-meson, 
however there is also a small contribution from 7t°7t ± 7t ± 
final states not related to 77 decays. Since the sum of the 
mass of the three pions sets a low mass limit and the 
total available energy a high mass limit of the distribu- 
tion, the background populates similar invariant masses 
as the 77 decays and cannot be efficiently suppressed by 
a cut on invariant mass. For further analysis all entries 
in the invariant mass spectrum have been accepted, but 
the resulting cross section has been corrected for the ob- 
vious background contribution at large invariant masses. 
For this correction the left hand side of the peak was mir- 
rored around the peak position (dashed histogram) and 
the access on the right hand side above this mirror image 
was assigned to the background. This procedure results 
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most likely in a slight underestimation of the total back- 
ground. 



4 Determination of cross sections 

The absolute normalization of the cross sections was ob- 
tained (see ref. [34] for details) from the target density 
((1.52±0.01) g/cm 2 ), the photon flux, the decay branch- 
ing ratio of tt° — > 77, and the acceptance and detection 
efficiency of the TAPS detector. The latter was obtained 
from MC simulations with the GEANT3 package [46] . For 
the simulations the pion pairs have been generated accord- 
ing to phase space distributions in quasifree kinematics 
taking into account the Fermi motion of the nucleons. FSI 
was taken into account as in [41) in a simple approxima- 
tion, where pions scatter randomly off nucleons according 
to their mean-free path. It was then checked, that missing 
mass-, kinetic energy-, and angular distributions were in 
agreement with the simulations. 

For the extraction of the rj photoproduction cross sec- 
tion from the r\ — > 77 decay the detection efficiency was 
obtained from MC simulations as function of laboratory 
angle and energy of the mesons as in [47] . This means that 
it becomes independent on assumptions about energy or 
angular distributions of the mesons. This is not possible 
for the 77 — > 7r°7r + 7r~ channel since the analysis used cuts 
on the reaction kinematic which depend on the incident 
photon energy. However, in this case the cross sections ex- 
tracted from the 77 — > 77 channel can be used as input for 
the event generator. This is possible since the relatively 
long-lived 77 mesons decay always outside the nucleus so 
that the decay pions are not influenced by additional FSI 
effects. The results for the total cross section of 77 produc- 
tion are summarized in fig. [8] The agreement between the 
previous data from [47] and the present results is quite 
good for both decay channels and can be used to estimate 
the systematic uncertainties. 



5 Systematic uncertainties 

The systematic uncertainties for the two reaction channels 
and for different energy ranges arise from different sources, 
which are briefly discussed in this section. 

We start with double ir° photoproduction below the 77 
production threshold (i.e. for E~ t < 550 MeV). The identi- 
fication of photons in the TAPS detector is very clean (see 
e.g. ref. [50]), and the suppression of background by the 
invariant mass analysis (see fig. [3]) is efficient. As shown 
in fig. [4] (left hand side), no background is visible in the 
missing mass spectrum for incident photon energies below 
500 MeV, and the shape of the missing mass spectrum 
is excellently reproduced by the Monte Carlo simulation. 
The missing mass cut is therefore irrelevant in this energy 
region. Furthermore, in a previous experiment [9], with an 
earlier stage of the TAPS detector, the reaction was inves- 
tigated below the 77 threshold from an analysis of events 
with three photons. In this case, only identification of one 
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Fig. 8. Total cross sections for r\ photoproduction off calcium. 
Compared are the results for the r\ — > 77 and r\ — > 7r°7r + 7r~ 
decays to the results from a previous measurement [47] . 



7T° meson and an additional photon was required. Such an 
analysis was possible since the only physical background in 
this energy range comes from the reaction 7 A — > j'it'^N. 
Compared to double n° production, the cross section of 
this background reaction is smaller by more than an or- 
der of magnitude [51], and the contribution to the three- 
photon sample was suppressed by another order of magni- 
tude relative to the 7r°7r° channel via detection efficiency 
(three-out-of-three versus three-out-of-four detected pho- 
tons). The result of this analysis for the total cross section 
is compared in fig. [9] to the present result. The systematic 
difference between the two excitation functions is only 3 
%, although the detector geometry, the cuts for the re- 
action identification, and the detection efficiencies were 
very different. We therefore estimate the systematic un- 
certainty of the double ir° data below the 77-threshold with 
5 %. 

At higher incident photon energies two effects increase 
the systematic uncertainty for the double ir° channel. Back- 
ground from r\ — ► 37r° decays, where the decay photons 
from one pion escape detection, becomes important. At 
the same time, the probability of re-scattering of the pions 
from the double 7r° production in the nucleus increases, 
since the mean free path of pions decreases strongly with 
increasing kinetic energy [9]. Together with the smearing 
of the reaction kinematics due to Fermi motion, the re- 
scattering results in a significant fraction of events which 
in the missing mass analysis overlap with the 77 back- 
ground (see fig. [4]). In the analysis, only events with miss- 
ing masses between 100 and -100 MeV were accepted for 
double 7T° production. The fraction of true 2ir a events, 
which do not pass this cut, was determined from the Monte 
Carlo simulations and corrected. The background from 77 
events which pass this cut can be subtracted, since the 77 
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Fig. 9. Total cross sections for double 7r° production extracted 
from events with three [5] and four detected photons. 



cross section is precisely known from the r\ — > 77 decay 
channel. This, however, involves a Monte Carlo simulation 
of r\ — > 37T° events. Both simulations contribute to the sys- 
tematic uncertainty. Among these two, the 2ir° simulation 
is more critical since it involves model assumptions about 
the pion FSI. From a variation of the missing mass cut we 
estimate the uncertainty at the 10 % level for the high- 
est incident photon energies. Therefore a linearly rising 
systematic error (0 % at 550 MeV, 10 % at 800 MeV) is 
added in quadrature. 

The systematic uncertainty for the mixed charge chan- 
nel is larger since for the charged pion we do not have the 
efficient invariant mass filter. A check of the simulated 
detection efficiency for charged pions is obtained from 
the analysis of the r) — > 7r°7r + 7r~ decay, where it enters 
squared. The agreement between the cross sections from 
the two different 77 decays (see fig. |HJ limits the systematic 
uncertainty for detection of a tt ^ pair to the 15 % range. 
At low incident photon energies the main source for sys- 
tematic uncertainties is background due to misidentifica- 
tion of recoil protons from single ir° production as charged 
pions (see sec. SJ, which is suppressed by a missing mass 
cut (see fig. [6]). Also from the comparison of the simulated 
missing mass spectrum to the measured distribution, a 
15 % uncertainty at the lowest incident photon energies 
appears realistic. It should be noted that the fraction of 
true double pion events, which do not pass the 150 MeV 
cut, enters of course via the Monte Carlo simulations into 
the detection efficiency Systematic effects therefore arise 
only from the correct representation of the double pion 
reaction in the simulation and from residual proton back- 
ground, which is not eliminated by the cut. The relative 
importance of the proton background decreases with in- 
creasing incident photon energy since the cross section ra- 
tio of single it to production drops by more than an 
order of magnitude between 450 and 600 MeV, although at 
the same time the kinematic separation of the background 
gets worse due to pion re-scattering. 

The relative contribution of background from ?7-decays 
(77 — > 7T°7r + 7r _ ) is smaller by roughly an order of magni- 



tude compared to the neutral channel. This is a combined 
effect of cross section ratios, r\ decay branching ratios and 
misidentification probability (ir° and one of the charged 
pions must be detected while any pair out of three tt° 's is 
possible) . Therefore it was not attempted to suppress this 
background with further cuts, but the simulated contri- 
bution was subtracted from the cross section. The related 
systematic uncertainty rises from threshold to high inci- 
dent photon energies to a few per cent and counteracts 
the dropping of the uncertainty related to the above dis- 
cussed background from single tt° production. In addition, 
at higher incident photon energies a small contamination 
from 7r°7r ± 7r ± events not originating from r] decays may 
contribute (few per cent level, see fig. [7J). We have thus 
assumed a conservative total uncertainty of 15 % for all 
incident photon energies. 

Finally, an overall normalization uncertainty for the 
luminosity (target surface thickness and photon flux) of 2 
- 3 % must be added. 



6 Results 

The total cross sections for the two double pion production 
channels are compared in fig. [TO] to data off the nucleon. 
Two remarks must be made. The unpolarized double pion 
data off the nucleon measured with DAPHNE at MAMI 
[2"Tll2"8"] was recently revised [32] . The data shown for the 
11(7, 7r°7r _ )p reaction are the result of this re-analysis [29] . 
The total cross section for production off 40 Ca from 

the TAPS experiment reported in [9] is slightly higher 
than the present result (approximately 10 % at 800 MeV). 
This discrepancy is due to the background from protons 
leaking into the pion TOF-versus-E band, which was not 
efficiently enough suppressed in the earlier analysis. 

When properly scaled by the mass number dependence 
(as discussed in detail in [9"l[TO]1 the total cross section 
for the 7T 7r° channel is in excellent agreement with the 
deuteron data over the full energy range, indicating that 
only trivial effects like FSI influence the nuclear cross sec- 
tion. The situation is less clear for the reaction. 
The nuclear data show a smooth behavior very similar 
to the n(7,7r°7r~)p reaction, while the p(~f,Tr°ir + )n exci- 
tation function has a pronounced peak in the second res- 
onance region. This behavior might contribute to the still 
not understood complete suppression of the second reso- 
nance bump in total photoabsorption data [52]. At least 
for this channel, a part of the effect might be due to a 
difference in proton and neutron cross section. However, 
since the 7r°7r + cross section is measured off the free proton 
and the 7r°7r _ cross section off the neutron bound in the 
deuteron a more detailed investigation of possible nuclear 
effects in the deuteron is highly desirable. 

The main motivation for this experiment was the in- 
vestigation of the 7r — 7r invariant mass distributions in 
view of possible in-medium modifications of the strength 
in the scalar, iso-scalar channel. For two reasons low in- 
cident photon energies are advantageous for this analysis. 
First this avoids the additional background from triple 



8 



F. Bloch et al.: Double pion photoproduction off Ca 



50 



40 



.O 30 



< 

D 20 



10 



■ 40 Ca(y,7i°7i +/ ")X 

a p(y,7i°7i + )n 

o n(Y,7tV)p 

• 40 Ca(Y,rcV)X 

o d(y,7i 7i°)np 



* .4 



<H 



«) «)•(« :» 



300 



400 



500 



600 



700 



800 



E y [MeV] 



Fig. 10. Total cross sections for and 7r°7r photopro- 

duction off 40 Ca compared to d(-y,-K ir°)np [36] respectively to 
p(7, 7r°7r + )n [37] and 71(7, 7r°7r~)p [281129] . Cross sections are 
normalized to A e ff = A 2 ^ 3 for calcium and A e ff = A for 
A = 1,2. 



MeV) are summarized in fig. [IT] and compared to previ- 
ous results for carbon [UJ. The distributions for the two 
different charge channels and also for the two different nu- 
clei are very similar. For 40 Ca we do not observe such an 
important shift of the 7r°7r° invariant mass distributions 
towards small invariant masses as has been found for lead 
[41] . neither relative to the ■k q tt ± channel for Ca nor in 
the comparison to 7r°7r° off carbon nuclei. 

For a more detailed comparison the cross section ratios 
C^{Ca/C) and C^ttV/ttV^ defined by: 



(T Ca (Tr ir )dM 



a c {Tr TT°)dM 
1 dacain ^) 



a C a(^ ^°)dM I £7 ao (7T 7T ± )dM 



(3) 
<4) 



where ax(^ ^ c ), V=carbon, calcium are the total cross 
sections and da/dM the invariant mass distributions, are 
shown in fig. 1121 In line with the previous results for car- 
bon and lead [41] some low mass enhancement of the 7r°7r° 
final state versus the 7r°7r ± final state is also visible for cal- 
cium. However, there is no significant difference between 
the 7r°7r° distributions off carbon and off calcium. This 
means that the size of possible in-medium effects is compa- 
rable for carbon and calcium, while significantly stronger 
effects had been found for lead |4"T1. 



pion decays of the r\ meson, which would give rise to addi- 
tional systematic uncertainties. Furthermore, at low inci- 
dent photon energies the kinetic energies of the produced 
pions are small. This is advantageous because low energy 
pions, which have only a small chance to excite the nu- 
cleoli to the A resonance, are much less affected by FSI 
[9]. On the other hand, the cross section drops strongly 
towards the production thresholds, so that some compro- 
mise is necessary. 
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■K°n ± are compared to each other and to the distributions 
from 12 C gTl. 
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Fig. 12. Cross section ratios C-n-niCajC) and 
Ctttt(ti Ti / n ti^) for low invariant masses. 



7 Comparison to results of the BUU-model 

The quantitative discussion of the experimental results re- 
quires a detailed understanding of the 'trivial' in-medium 
effects like the smearing of cross sections due to the mo- 
mentum distributions of the bound nucleons, collisional 
broadening of nuclcon resonances due to additional decay 
channels like NN* — > NN, Pauli-blocking of final states, 
and final state interactions of the pions. These effects 
have been intensively studied in the framework of cou- 
pled channel transport models based on the semi-classical 
Boltzmann-Uehling-Uhlenbeck (BUU) equation (see [53] 
54,55] for details of the model). Recently Miihlich et al. 
[5T>] and Buss et al. [S7] have performed detailed calcula- 
tions of the nuclear double pion photoproduction reactions 
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Fig. 13. Total cross sections compared to results of the BUU 
model |57j . The bars at the bottom represent the systematic 
uncertainty of the data (see sec. [5}, the dashed lines represent 
the error band for the BUU calculation. 
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Fig. 14. Pion-pion invariant mass distributions compared to 
results of the BUU model |57j . The bars at the bottom rep- 
resent the systematic uncertainty of the data (see sec. [5}, the 
dashed lines represent the error band for the BUU calculation. 



in the framework of this model. In the latter work, special 
emphasis was laid on the description of the scattering and 
absorption properties of low energy pions. 

Apart from the inherent assumptions and approxima- 
tions made in semi-classical transport models, systematic 
uncertainties of the calculations arise from the uncertainty 
of the cross sections of the elementary production reac- 
tions and the uncertainty of the cross sections involved in 
the description of pion FSI. Some indication for the cor- 
rect treatment of pion FSI comes from a comparison [58j 
of model calculations and experimental results for pion- 
induced double charge exchange reactions, which are very 
sensitive to details of pion FSI. It is thus assumed, that the 
major uncertainty is related to the experimentally deter- 
mined elementary production cross sections of pion pairs 
off nucleons. Buss and co-workers |57j estimate a system- 
atic 20 % uncertainty of their results for the 7r°7r° channel 
from a simultaneous variation of all input cross sections 
to their ±ler limits. This estimate is conservative since no 
cancellation of systematic errors from different reaction 
channels is considered. The largest uncertainty is related 
to the 771 — > nir°ir° channel. The uncertainty estimated 
in [57] for the ir ^ final state is much smaller, however 
this was based only on the statistical uncertainties of the 
input data and has neglected the systematic ones, which 
are also typically on the 10 - 15 % level. Therefore we as- 
sume also for this channel a 15 % systematic uncertainty. 
These model uncertainties and also the systematic uncer- 
tainties of the data are indicated in figs. 1131141 as dashed 
lines, respectively bar charts. 

The overall agreement between the model results and 
the data for the total cross sections is good within the 
systematic uncertainties (see fig. [T3")) , it is excellent for the 
7r°7r° channel in the energy range of interest for the invari- 
ant mass distributions. Up to now, reliable BUU calcula- 
tions are not available for incident photon energies above 
700 MeV, since in this region the elementary free nucleon 



cross sections have too large uncertainties for some chan- 
nels. 

An important result of this model calculations is, that 
even without any explicit in-medium modification of the 
7r°7r° channel, the respective invariant mass distributions 
show a significant softening for heavy nuclei [57J. In the 
model, this effect arises solely from FSI of the pions. Part 
of it comes from the fact that the irN absorption cross sec- 
tion increases strongly with pion kinetic energy, so that pi- 
ons with larger kinetic energies are more strongly depleted 
via the ttN — > A, AN — > NN reaction path. But not only 
the absorptive part of the FSI is important, also impor- 
tant are re-scattering processes which tend to decrease 
the pion kinetic energy and thus the average pion - pion 
invariant masses. This effect is enhanced due to charge 
exchange scattering, which mixes the contributions from 
the different charge channels. Since the total cross section 
for 7r°7r production is much larger than the 7r°7r° cross 
section, the latter receives significant side feeding from the 
mixed charge channel via tt^N — > AV scattering, which 
increases the fraction of re-scattered low-energy pions in 
this channel. In the same way, re-scattering of ir + Tr~ con- 
tributes to the tt ^ channel. This result means, that the 
relative softening of the ir°ir° invariant mass by itself is 
no evidence for an in-medium modification of the pion - 
pion interaction, but the argument becomes quantitative, 
requiring a comparison of the experimental results to cal- 
culations which account for the trivial FSI effects. 

The invariant mass distributions for the two energy 
ranges are compared to the BUU results in fig. [14] The 
general tendency of a softening of the distributions is re- 
produced by the calculations as a consequence of FSI ef- 
fects. On a quantitative level, the data for the ir°ir° invari- 
ant masses for the lowest incident photon energies seem 
to be slightly stronger downward shifted than predicted. 
One should note that this effect cannot be completely ex- 
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Fig. 15. Cross section ratio C 7 r7r(vr 7r /7i"°7r ± ) compared to 
the results of the BUU model [57]. Symbols: data, curves BUU 
results [57]. Left hand side: incident photon energies 400 - 500 
MeV, right hand side: incident photon energies 500 - 550 MeV. 



plained by the systematic uncertainties of data and model 
results, since the uncertainties are not strongly invariant 
mass dependent but rather related to the absolute scale of 
the cross section. The cross section for the mixed charge 
channel is somewhat overestimated by the model, in par- 
ticular for large invariant masses. This absolute normaliza- 
tion cancels in the C 7 r T (7r 7r /7r 7r ± ) ratio, which is com- 
pared to the BUU results in fig [15] Within the statistical 
uncertainties, the agreement between data and model re- 
sults is good. This means that the relative shapes of the 
7r°7r° and ir ^ invariant mass distributions can be almost 
completely reproduced from FSI effects. 



8 Summary and conclusions 

Total cross sections and pion-pion invariant mass distri- 
butions have been measured for 7r°7r° and -k ^ photo- 
production off 40 Ca nuclei. When scaled by A 2 / 3 the total 
cross sections agree quite well with the average of the el- 
ementary cross sections off the proton and the neutron. 
The only exception is the mixed charge final state in the 
second resonance region, where the reaction off the free 
proton shows a more pronounced signal for the resonance 
bump. 

The invariant mass spectra show a similar effect as 
already reported in [41] for carbon and lead nuclei, namely 
a softening of the 7r°7r° distributions relative to the ir^ir^ 
distributions. The strength of the effect is comparable to 
carbon. 

The data have been compared to calculations in the 
framework of the BUU model [57] . A sizable part of the 
in-medium effects can be explained by the model by final 
state interaction effects, which tend to shift re-scattered 
pions to smaller kinetic energies. Only for the lowest inci- 
dent photon energies a small additional downward shift of 
the strength to small invariant masses for the 7r°7r° chan- 
nel may be visible. 
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